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ABSTRACT 

We show that in a supernova core the annihilation process v c D c — > v^ T v^ r is always more 
important than the traditional reaction e + e~ — ► v^ T v^ r as a source for muon and tau neutrino 
pairs. We study the impact of the new process by means of a Monte Carlo transport code with 
a static stellar background model and by means of a self-consistent hydrodynamical simulation 
with Boltzmann neutrino transport. Nucleon bremsstrahlung NN — > NNz^. r £y T is also included 
as another important source term. Taking into account v c v c — > v^.tV^.t increases the and v T 
luminosities by as much as 20% while the spectra remain almost unaffected. In our hydrodynam- 
ical simulation the shock was somewhat weakened. Elastic v^^v^ and v^. T Ve scattering is not 
negligible but less important than v^e scattering. Its influence on the v^^ fluxes and spectra 
is small after all other processes have been included. 

Subject headings: neutrinos — supernovae: general 



1. INTRODUCTION 

The treatment of and v r transport in numer- 
ical supernova (SN) simulations has been some- 
what schematic in the past. However, with the 
advent of numerical Boltzmann solvers for the neu- 
trino transport (Mezzacappa & Bruenn 1993, Mez- 
zacappa & Messer 1999; Yamada, Janka, & Suzuki 
1999; Burrows et al. 2000; Rampp & Janka 2002) 
and their application to the post-bounce phase 
of stellar core-collapse models (Rampp & Janka 
2000; Mezzacappa et al. 2001; Liebendorfer et al. 
2001) a new level of accuracy has been achieved. 
Evidently it is desirable that in consistent state- 
of-the-art simulations the uncertainties are not 
dominated by overly crude approximations of the 
microphysics which governs the neutrino interac- 
tions. For example, it has been recognized that nu- 
clear many-body correlations (Burrows & Sawyer 
1998; Reddy et al. 1999) or weak-magnetism ef- 
fects in neutrino-nucleon interactions (Vogel & 
Beacom 1999; Horowitz & Li 2000; Horowitz 2002) 



should be included. 

Previous simulations used iso-energetic scatter- 
ing on nucleons i^N — > Nf M as the main opacity 
source for transport, elastic scattering on elec- 
trons and positrons v^e^- — > e^v^ as the main 
energy-exchange reaction, and e + e~ — > v^v^ as 
the only source term for production. (Here and 
in what follows we use symbolically for either 

or i/ T .) However, it is now generally accepted 
that nucleon bremsstrahlung NN — » NNz^P^ is 
important or even dominant as a neutrino source 
reaction (Suzuki 1991,1993; Hannestad & Raffelt 
1998; Thompson, Burrows, & Horvath 2000), and 
that nucleon recoils have a significant impact on 
the emerging flux spectrum (Janka et al. 1996; 
Raffelt 2001; Keil, Raffelt, & Janka 2002). 

In this paper we show that in addition v e v e — > 
v^v^ and its inverse reaction should be included 
because this process is always far more important 



than e + e 



as a neutrino source (Fig. 1). 



Conversely, v^v c and v^v c scattering turns out to 
be less important than v^e^- scattering and thus 
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Fig. 1. — Feynman graph for the annihilation pro- 
cesses producing v^v^ pairs. 

not crucial for the flux and spectra formation. 
Bond (1978) has previously discussed these pro- 
cesses in the context of neutrino transport in su- 
pernovae, but no assessment of their relative im- 
portance has been performed. 

It is straightforward to see that v e v e — > v^v^ 
must be important in a SN core. The region of 
interest is well inside the transport sphere, i.e. 
interior to the radius where v^N — > Ni^ scatter- 
ing freezes out. This region is always deeper in 
the star than the freeze-out radius of the charged- 
current reactions v c n «-> pc~ and 9 c p <-> ne + , 
implying that in the trapping region v c and 
v c are essentially in local thermodynamic equi- 
librium (LTE). Therefore, in this region the pro- 
cess v c v c <-» v^Ufj, will serve to create or destroy 
i^n^n pairs in the same way as the traditional 
e + c~ <-> v^v^ process. With vanishing chemical 
potentials of c~ and v e , the rate of the new process 
turns out to be about twice that of the traditional 
one so that the combined source strength would 
be three times the traditional one. 

In the presence of non-vanishing chemical po- 
tentials the enhancement is even larger. The pair 
production rate due to e + e~ — > v^v^ is shown in 
Fig. 2 as a function of the degeneracy parameter 
r] e = Hc/T for the electrons, and likewise, the rate 
for v c v c — ► VpVp as a function of t\ Vb . The produc- 
tion rate is a decreasing function of rj. Since in 
the relevant regions of the SN core r\ Va < r/ c , the 
traditional process of e + c~ annihilation is reduced 
more strongly than the new v c v c rate. The latter 
therefore dominates even more. 

In order to estimate the impact of the new 
source reaction on the neutrino fluxes and spectra 
we conduct two separate numerical investigations. 
First, we perform a Monte Carlo transport simu- 
lation on a static stellar background model. While 
this approach has the disadvantage of not follow- 



ing the evolution of the neutron star atmosphere 
self-consistently, it enables us to disentangle the 
individual effects of different neutrino processes on 
the transport and spectra formation in a system- 
atic study. Next, we perform two full-scale numer- 
ical post-bounce simulations where once the new 
process is included and once it is left out. This al- 
lows us to verify that the effects established by the 
Monte Carlo results are generic and also show up 
in self-consistent radiation-hydrodynamical mod- 
els. The use of these two different approaches and 
independent codes also helps making sure that our 
results do not depend on details of the technical 
implementation or the particular numerical reso- 
lution. 

We begin in Sec. 2 by comparing the e + e~ and 
v c v c reactions. In Sec. 3 we discuss the results of a 
Monte Carlo study of neutrino transport while in 
Sec. 4 we describe the self-consistent hydrodynam- 
ical simulations coupled with a Boltzmann trans- 
port solver. We summarize our findings in Sec. 5. 

2. ELECTRON VS. NEUTRINO PAIR 
ANNIHILATION 

We begin by comparing the two pair annihila- 
tion processes 

^o + ^o -» Vp + Vp, (1) 

e+ + e~ — ► Vn + Dp. (2) 




Fig. 2. — Pair production rates by the process 
v c v c —> v^v^ as a function of rj Ve (upper line) and 
e + e~ — * VpVn as a function of rj e (lower line). We 
used T = 12 MeV and ^ = 0. 
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Table 1 
Weak interaction constants. 



Process: 


0+C <-> V^Vp 


Z/ e I/ e <- 




Cv 


-\ + 2 sin 2 6>„, 






C A 


i 

2 







In the relativistic limit where the electron mass 
can be neglected, their squared and spin-summed 
matrix elements are of the form 

]T|-M| 2 = 8 G 2 F [ (Cy + Ca) 2 u 2 

spins 

+ (C v - C A ) 2 t 2 } , (3) 

with the Mandelstam variables t = — 2k\ ■ k 3 and 
u = —2ki ■ fc 4 . The momenta are assigned to the 
particles as indicated in Fig. 1. The matrix ele- 
ments for the two processes differ only in the weak 
coupling constants Cv,a shown in Table 1. 

The pair production rate is obtained by ap- 
propriate phase-space integrations, including par- 
ticle distributions and blocking factors (Yueh & 
Buchler 1976; Hanncstad & Madsen 1995). We 
use n Vfl = even though there could be a small 

chemical potential due to a non-vanishing con- 
centration of muons in the core and due to dif- 
ferent transport properties of and v^. As- 
suming that c ± , v e , and D c are all in LTE, the 
rates of pair creation vs. degeneracy parameters 
r] c and ?7„ o are shown in Fig. 2. In the dense re- 
gions of a SN core below the neutrino spheres the 
phase space distribution of electron neutrinos is a 
Fermi-Dirac function with degeneracy parameter 
f]v e = Vc + Vp ^ Vn < Vc so that the new process is 
always more important than e + e~ annihilation. 

An interesting difference between the two pro- 
cesses arises in the differential production rates, 
i.e. the and production rates as functions of 
neutrino energy e. As a first case we take 77 = 
for both e and v c and show the differential produc- 
tion rate d 2 n/dedt in Fig. 3. For both processes 
the differential rate is the same for and v^, i.e. 
they are produced with the same spectra. 

However, in general the "parent particles" will 
have a significant chemical potential. Taking r] e = 



r] Ve = 10 we show the differential production rates 
in Fig. 4. In case of the e + e~ process (upper 
panel) the differential rates are similar for and 
Vp. This is understood by the fact that the values 
of (C v + CO 2 =s 0.54 2 and (C v - C A ) 2 0.46 2 
are quite similar so that the u 2 and t 2 terms in 
Eq. (3) are almost equally important. Therefore, 
interchanging and P M , corresponding to an ex- 
change of u and t, has no big effect. 

This is not true for the neutrino reaction, where 
(C V + C A ) 2 = 1 and (C V -C A ) 2 = and thus only 
u 2 contributes. Replacing u by t now changes the 
kinematics of the process, which in turn modifies 
the rate if the distribution of v c differs from that of 
v e , i.e. if r) Ve ^ 0. The and P M spectra shown in 
the bottom panel of Fig. 4 are indeed very differ- 
ent from each other, although the total production 
rates of and are, of course, equal. 

One can easily understand why on average 
have larger energies than v^. We first look at 
v c v c — » v^Vfj, in the center of momentum (CM) 
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Fig. 3. — Differential and production rates 
d 2 n/dedt vs. neutrino energy for i] c = r\ Ve = and 
T = 12 MeV. 



3 




20 40 60 80 100 120 140 160 180 
Energy [MeV] 

Fig. 4. Differential and P M production rates 
d 2 n/dedi for rj e = r) v<s = 10 and T = 12 MeV. 
Upper panel for e + c~ — > t^f^, lower panel for 

frame. The differential cross section is 
dcr 



, „ ^i e 2 (l + C os0) 2 , (4) 
dcos6» 47r 

where 6> is the angle between the ingoing v c and the 
outgoing f^, or equivalently, between the ingoing 
v e and the outgoing v^. Put another way, forward 
scattering is favored and backward scattering for- 
bidden. This is due to angular momentum conser- 
vation. The ingoing v c and v e have opposite helic- 
ities and, in the CM frame, opposite momenta, so 
that their combined spins add up to 1. The same 
is true for the outgoing particles so that backward 
scattering would violate angular momentum con- 
servation. In the rest frame of the medium the 
ingoing v c tends to have energies of the order of 
its Fermi energy, while the ingoing P e tends to have 
energies of order T. Because forward scattering is 
favored, the outgoing tends to inherit the larger 
energy of the ingoing v e . 

The differences of the source spectra, however, 



do not translate into significant spectral differ- 
ences of the Ufj, and fluxes emitted from the 
SN core. While pair annihilations and nucleon 
bremsstrahlung are responsible for producing or 
absorbing neutrino pairs and thus their equilibra- 
tion with the stellar medium below the "neutrino- 
energy sphere," other processes, notably v^e* 1 
scattering and nucleon recoils, are more efficient 
for the exchange of energy between neutrinos and 
the medium between the equilibration and trans- 
port spheres. In our numerical runs we will find 
in fact that adding the new process to a SN simu- 
lation primarily modifies the flux with only minor 
modifications of the spectrum. 

If v c v c — > VpVp is important relative to e + e~ — > 
Vypy, one may wonder if processes of the form 



(5) 
(6) 



could be of comparable importance as f M e± scat- 
tering. Figure 5 shows the rates for f M scattering 
on v e and 9 e , and those for scattering on e+ and 
e~ as functions of r} v<s and r] e , respectively. The 
rates are normalized to the rate at i] c = 0. 

In contrast to the annihilation rates, the scattering 
rates rise monotonically with r\. Therefore, even 
though neutrino-neutrino scattering would domi- 
nate if all chemical potentials were zero, for realis- 
tic situations with rf Ve < rj c we expect that scatter- 
ing on e"* 1 has 1-2 times the rate of scattering on 
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Fig. 5. — Thermally averaged scattering rate for 
Ujx on c^ as a function of rj c (lower line) and for 
on i/ e and v c (upper line) as a function of r\ Vi! . 
The rates are normalized to the scattering rate on 
e ± at r} c = 0. We used T =12 MeV and r\ v = 0. 
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v e and v e . Therefore, neutrino- neutrino scattering 
is expected to be a relatively minor correction. In 
our Monte Carlo studies we will indeed find that 
this process has only a small effect on the neutrino 
spectra and fluxes. 

3. MONTE CARLO STUDY 

To study the impact of the new annihilation 
process on the fluxes and spectra we first use a 
Monte Carlo method for neutrino transport on a 
static background stellar model. For this purpose 
we have adapted the Monte Carlo code of Janka 
& Hillebrandt (1989a, b) to include additional pro- 
cesses and nucleon recoil. We use this code for 
an extensive parameter study of spectra for- 
mation that will be documented elsewhere (Kcil 
ct. al. 2002). 

As a stellar background we employ a model 
originally provided to us by B. Messer for an ear- 
lier study of neutrino spectra formation (Raffclt 
2001). Using this model again for our present 
study facilitates a comparison with this previous 
work. The model is based on a full-scale Newto- 
nian collapse simulation of the Woosley & Weaver 
15 Mq progenitor model labeled sl5s7b. We use a 
snapshot at 324 ms after bounce when the shock 
wave is at a radius of about 120 km, i.e. the SN 
core still accretes matter. In Fig. 6 the tempera- 
ture profile is represented by the steps in terms of 
(e) ~ 3.15 T for each radial zone, i.e. by the aver- 
age energy of nondegenerate neutrinos (rj Vfl = 0) 
in LTE with the stellar medium. 

In Fig. 6 we also show the thcrmalization depth 
-Rtherm for several processes as a function of neu- 
trino energy e. The formal definition of i?thcrm 
in terms of an effective mean free path for energy 
exchange is given by the condition (cf. Shapiro & 
Teukolsky 1983; Suzuki 1989; Raffelt 2001) 




(7) 



for the effective optical depth r e s of a partic- 
ular equilibrating process with mean free path 
Aj (e) among the opacity producing reactions hav- 
ing mean free paths Xj(e). The solid line in Fig. 6 
is for v e n — ► pe~, i.e. it represents the energy- 
dependent v e sphere. The dashed line is the analo- 
gous P e sphere due to D c p — > ne+ . Finally, the dot- 
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Fig. 6. — Thermalization depth -Rthcrm (hori- 
zontal axis) for different processes as a function 
of neutrino energy e (vertical axis) in our back- 
ground model. The energy-exchanging processes 
are v c n — ► pc~ (solid line), P e p — > ne + (dashed 
line), and v^v^ annihilation to v e 9 e pairs (dotted 
line). The steps represent the temperature profile 
of the stellar model in terms of (e) ~ 3.15 T. 



ted line shows where or of given energy last 
participate in the v^v^ — > v e D e process, assuming 
that the annihilation partners are distributed ac- 
cording to LTE. Put another way, the dotted line 
is the energy-dependent freeze-out sphere for our 
new process. It is always at much smaller radii 
than the v c and D c spheres. Therefore, in those 
regions where VuPu, *-> v^v^ is effective we may as- 
sume LTE for v c and v c . 

We therefore implemented the new process in 
our Monte Carlo code by using the same subrou- 
tine as for the e + e~ pair process, except for in- 
serting the appropriate weak interaction constants 
Cy,A and replacing r] e by r\ Ve . 

We summarize the characteristics of the emerg- 
ing neutrino flux for our runs in Table 2. As intro- 
duced by Janka & Hillebrandt (1989a) we charac- 
terize the neutrino spectrum by its first two energy 
moments (e)fl ux and (e 2 )fl ux and define the "pinch- 
ing parameter" by following Raffclt (2001) as 



Pflu 



1 (e 2 )fl, 



a <e} 2 



flux 



with the energy moments 

J oo de/+ 1 1 d M / !y (e, A1 )e"+ 2 M 



(Ofl. 



/o TOde /-i d M/«>(e,M)e 2 M 



(8) 



(9) 
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Table 2 

Spectral characteristics of neutrino fluxes from Monte Carlo transport. 



Energy exchange 


(e)flux 


(e 2 )flux 


Pflux 


T 


V 








[MeV] 


[MeV 2 ] 




[MeV] 




[10 51 


transport 
















original run 




17.5 


388. 


0.97 


5.2 


1.1 


14.4 


s p 




16.6 


362. 


1.01 


5.3 


-0.3 


15.8 


s p 


n 


16.9 


369. 


0.99 


5.3 


0.4 


20.2 


b r s p 




14.2 


255. 


0.98 


4.2 


1.1 


14.8 


b r s p 


n 


14.4 


264. 


0.97 


4.3 


1.2 


17.6 


b r s 


n 


14.4 


263. 


0.97 


4.3 


1.2 


17.0 


b r sn p 


n 


14.3 


260. 


0.97 


4.3 


1.2 


17.9 


Dfj, transport 
















s p 


n 


16.9 


368. 


0.99 


5.2 


0.6 


20.6 


b r s p 


n 


14.4 


263. 


0.97 


4.2 


1.3 


17.8 


b r s 


n 


14.4 


262. 


0.98 


4.3 


1.1 


16.8 



Note. — For energy exchange, "b" refers to bremsstrahlung, "r" to recoil, 
"s" to scattering on electrons and positrons, "p" to e + e~ pair annihilation, 
"n" to v a v e pair annihilation, and "sn" to scattering on both, and u e , v c . 
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of the emergent flux spectrum. Here f v (e,n) is 
the neutrino distribution function in energy-angle 
space with fi being the cosine of the angle of neu- 
trino propagation relative to the radial direction. 
The constant a for a Fermi-Dirac distribution at 
zero chemical potential is 



a 



• 486000 C3C5 



49 7T 8 



1.3029. 



(10) 



Then p = 1 signifies that the spectrum is ther- 
mal up to its second moment, while p < 1 signifies 
a pinched spectrum (high-energy tail suppressed) , 
and p > 1 an anti-pinched spectrum (high-energy 
tail enhanced). For p < 1.023 it is common to ap- 
proximate the spectrum as a nominal Fermi-Dirac 
distribution characterized by a temperature T and 
a degeneracy parameter 7/ which are chosen such 
that (e)fl U x and (e 2 )fl ux are reproduced (Table 2). 
Finally we show the neutrino luminosity in the last 
column of Table 2. 

We always include elastic z^N scattering which 
provides the dominant opacity contribution. In 
a given run we additionally include those energy- 
exchanging processes which are indicated in the 
first columns of Table 2. We use "b" to in- 
dicate bremsstrahlung NN — ► NNz^zy^, "r" for 



recoil in i/^N 



Nza, 



; 's" for v^e scattering, 
"p" for the traditional pair annihilation process 
e + c~ — > v^Vp, "n" for the new neutrino annihila- 
tion process v c D c — > v^v^, and "sn" for v^e^ 1 plus 
v^Vf, and z/ M z/ c scattering. 

In the first row of Table 2 we show the origi- 
nal flux characteristics of our background model 
from the calculation of Messer et al. (2002). Run- 
ning our code with the same neutrino reactions — 
scattering on c ± (s) and e+c - pair annihilation 
(p) — we find the results in the second row. We 
attribute the small differences to the different nu- 
merical approaches. In particular, the Boltzmann 
solver of Messer et al. (2002) works in practice 
with a limited number of energy bins. Moreover, 
there may be differences in the implementation of 
the microphysical reactions. Finally, in the Monte 
Carlo code we assume that neutrinos are in ther- 
modynamic equilibrium with the stellar medium 
at the inner boundary. This choice of bound- 
ary condition may have a small effect on our re- 
sults since the thermalization depth of c+c - pair 
processes is strongly energy dependent and thus 
low energy are forced into equilibrium by our 



boundary condition instead of the pair process. 
We interpret the first two rows of Table 2 as agree- 
ing satisfactorily well with each other. 

Next we add u e D e annihilation (n). The spec- 
trum remains almost unchanged, but the luminos- 
ity increases by about 30%. Therefore, the new 
process has a rather significant impact on the pre- 
dicted Vfj_ luminosity. 

However, other processes are also important 
which in the past have not been included in nu- 
merical simulations. Therefore, we switch off the 
new process and instead include bremsstrahlung 
(b) and recoil (r). Compared to the original run 
we obtain almost the same luminosity, but signif- 
icantly lowered spectral energies. Now we again 
include v c v c annihilation (n) and find that the 
spectra remain unaffected, but the luminosity in- 
creases by about 20%. Therefore, even with all 
other energy-exchanging reactions included, the 
f e z> e process still has an important effect on the lu- 
minosity. We finally switch off the traditional pair 
process (p), but keep the new one. The spectra 
remain unaffected, the luminosity slightly drops. 
It is evident that the v e 9 e process is by far the 
dominant leptonic source reaction for muon neu- 
trinos. Its importance relative to bremsstrahlung 
will depend sensitively on the background model 
(Keil et. al. 2002). 

As a last step we include the scattering on v c 
and 9 e in addition to all other processes. The rate 
for this process is typically half as large as the rate 
for scattering on c^, in agreement with Fig. 5 if 
we use T) c ~ 3 and ~ 0.3. The effect on the 
Vfj, flux and spectrum is minimal and in fact below 
the numerical resolution of our Monte Carlo runs. 

In the second part of Table 2 we finally show 
several runs for the transport of anti-neutrinos. 
Recall that v c D c — > v^v^ generates different source 
spectra for and v^. Of course, the small differ- 
ences between z/ M and scattering off electrons 
and positrons are also taken into account, as well 
as the small differences in e + e~ pair annihilation. 
Comparing the z/ M runs with those for we find 
excellent agreement. Therefore, the detailed spec- 
tral distribution of the pair rate is not important, 
only the total rate of absorption and production 
°f v^v^ pairs matters. 
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4. HYDRO DYNAMICAL SIMULATION 
WITH BOLTZMANN TRANSPORT 

The large flux increase found by the Monte 
Carlo method may not persist in a self-consistent 
hydrodynamical treatment where the stellar model 
can adjust in response to the modified transport. 
For this reason we have performed a Boltzmann 
transport simulation coupled with a full hydro- 
dynamics code (Rampp & Janka 2002) for the 
15 Mq progenitor model sl5s7b2 (Woosley, per- 
sonal communication; Woosley & Weaver 1995). 

In order to minimize the required modifica- 
tions of the code we treat the transport of 
and P M identically. Therefore, we use an average 
source strength for the and production from 
v e v e annihilation, ignoring the spectral differences. 
The Monte Carlo results suggest that this approx- 
imation is well justified. 

For computing the interaction rates we again 
assume LTE for v e and P e in all regions where 
the new annihilation process is important. This 
assumption also justifies introducing an energy 
source term for the new process directly into the 
medium energy equation in perfect analogy to the 
source term of the process e + e~ <-> vv. Similarly, 
the scattering reactions of and off v c and 
D c can be treated in full correspondence to the 
scattering off electrons and positrons by a simple 
change of the weak coupling coefficients. Thus a 
direct coupling of the and v c sectors of the neu- 
trino transport code is avoided. For the new pro- 
cess this is achieved only indirectly via the stellar 
medium as an intermediary 

Our baseline for comparison is a Newtonian 
simulation which includes the transport of neu- 
trinos of all three flavors. Neutrino-medium inter- 
actions for all relevant processes are implemented, 
notably nucleon bremsstrahlung, as described by 
Rampp & Janka (2002). Note, however, that nu- 
cleon recoils are not taken into account in the re- 
sults shown for the baseline simulation (thin lines 
in Figs. 7, 8 and 9). Rather, the charged-current 
and neutral-current neutrino reactions with nucle- 
ons are handled in the "old standard approxima- 
tion" of infinitely massive nucleons at rest follow- 
ing the treatment by Brucnn (1985) and Mezza- 
cappa & Bruenn (1993). 

In a second simulation we included the new lep- 
tonic process for v^v^ pair production and anni- 



hilation. The results for this run are depicted as 
thick lines in the figures. In Figs. 7 and 8 we show 
the evolution of the neutrino luminosities and of 
the root mean squared (rms) energies as functions 
of time after bounce. The rms energy is defined in 
the usual way (e.g. Messer et al. 1998) by 



(4 



Jo de Li d v.u( e ^) e5 



(11) 



i.e. using the neutrino energy distribution as a 
weight function. The results are given for an ob- 
server who is comoving with the stellar fluid at a 
radial location of 500 km. The Dopplcr-blucshift 
due to the infall velocity of the matter is rather 
small at this distance from the neutrino emitting 
neutron star so that the quantities are close to the 
observable properties at infinity. 

As expected, the luminosity is increased by 
10-20%, but the spectrum remains essentially un- 
affected. The enhanced energy loss leads to a 
somewhat faster proto-neutron star contraction. 
This causes small changes also in the region where 
the electron neutrino and antineutrino fluxes are 
built up and where these neutrinos finally decouple 
from the stellar background. As a consequence, 
the mean spectral energies of v c and D c are sys- 
tematically higher by a small amount (Fig. 8). 
Initially, the v c luminosity is also slightly larger 
(Fig. 7) but after about 150 ms post bounce the 
luminosities of v c and P c drop below the level of 
those in the reference simulation. This is caused 
by the decrease of the radii of the neutrino spheres 
in response to the accelerated contraction of the 
proto-neutron star. At 200 ms after bounce the 
v^.r luminosities of the two models then become 
nearly equal, probably as a consequence of two 
competing effects which seem to essentially com- 
pensate each other at later times: On the one hand 
side the neutrino emission of the nascent neutron 
star decays faster with time when the new process 
is included (as visible from the v e and v e luminosi- 
ties) , on the other hand the new process raises the 
energy loss in v^, T compared to v c and D c . 

The shock positions in both simulations evolve 
identically until about 100 ms after bounce. Then 
the shock is somewhat weakened by introducing 
the new reaction and expands only to a maxi- 
mum radius of 230 km instead of 250 km (Fig. 9). 
This can be understood again by the more com- 
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pact proto-neutron star, which causes higher in- 
fall velocities in the region of heating by v e and 
D e absorption behind the shock. This reduces the 
integral energy deposition by neutrinos as well as 
the pressure behind the shock front. 




50 100 150 200 250 300 350 
Time [ms] 



Fig. 7. — Evolution of the luminosities for v c 
(solid), D e (dashed) and (dotted) as functions of 
time after bounce for two hydrodynamical simula- 
tions with three-flavor Boltzmann neutrino trans- 
port. Thick lines represent results of a model with 
the new process included, thin lines of a reference 
calculation without this process. The luminosities 
are measured by an observer comoving with the 
stellar fluid at 500 km. 



30 u- 




Time [ms] 

Fig. 8. — Evolution of the rms energies for v c 
(solid), v c (dashed) and (dotted) as measured 
by an observer comoving with the stellar fluid at 
500 km. 

Since wc completed this study we have imple- 
mented the neutrino annihilation process in simu- 
lations with an approximate treatment of general 




Time [ms] 

Fig. 9. — Evolution of the shock positions. 

relativity (Rampp & Janka 2002). In these new 
simulations neutrino-nucleon interactions are im- 
plemented in terms of dynamical structure func- 
tions for correlated nuclear matter and also in- 
clude weak magnetism corrections (Rampp et al. 
2002). Put another way, these simulations include 
nuclear recoil as well as nuclear correlation effects. 
The impact of the new process in these more com- 
plete simulations is comparable to what we have 
found here. 

5. CONCLUSIONS 

We have shown that v e P e annihilation is the 
dominant leptonic source for v^^v^.r pairs in a 
SN core, far more important than the traditional 
c + c~ pair annihilation. Its importance relative 
to the nucleon bremsstrahlung process, which also 
has not been included in previous SN simulations, 
depends on the stellar profile. In our Monte Carlo 
studies with a background model representing the 
SN core during the accretion phase, the new pro- 
cess enhanced the luminosity by about 20% in a 
calculation where nucleon bremsstrahlung and en- 
ergy exchange by nucleon recoils were taken into 
account. Without these other new effects, v c v c 
annihilation has an even larger impact. 

In self-consistent hydrodynamical simulations 
with Boltzmann neutrino transport we find that 
during the first 150 ms the effect is similar to that 
obtained in the static Monte Carlo simulations and 
has a noticeable influence on the stellar evolution 
and structure. Later, the and luminosities 
approach those of the model without the new pro- 



9 



cess of v e v c annihilation. At that time, however, 
the v e and v e luminosities are smaller. Throughout 
the simulation the new reaction therefore works in 
the direction of making the v c and v c luminosities 
more similar to those of and P M . In the model 
with the new process the shock is somewhat weak- 
ened and reaches only smaller radii. 

In all of our Monte Carlo and hydrodynamical 
runs we confirm the naive expectation that the ef- 
fect of the new pair production and annihilation 
process on the neutrino spectra is minimal. The 
crossed process of v^ iT scattering off v c and P e also 
turned out to have a negligible impact on the emit- 
ted spectra. 

We conclude that state-of-the-art SN simula- 
tions should include the v c D c annihilation reac- 
tion to f^i-ty^i- pairs. While the effects of this 
process on the neutrino luminosities and spectra 
and on the shock propagation are not dramatic, 
they are nevertheless noticeable and not negligible, 
even after nucleon bremsstrahlung and nucleon re- 
coils have been included. Implementing the new 
process is not more difficult and not more CPU- 
expensive than the traditional e + e~ process. 
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